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I.  MANAGEMENT  REPORT 


A.  Summary 

The  major  development  in  the  last  few  months  has  been  the  demonstration 
of  broadband,  Sezawa  mode  ZnO  devices.  We  have  demonstrated  devices 
with  a  40  MHz  bandwidth,  i.e.,  approximately  25%  .  We  have  also  developed 
an  entirely  new  program  for  the  design  of  matching  networks  to  arbitrary 
loads.  We  believe  that  this  is  extremely  important  because  it  is  not 
based  on  the  usual  matching  network  approach.  Instead,  it  is  based  on 
the  development  of  our  adaptive  filter  concepts  in  which  we  postulate  the 
form  of  the  matching  network,  then  adjust  its  parameters  on  the  computer 
in  a  manner  equivalent  to  trimming  the  elements  to  obtain  the  optimum 
broadband  match.  The  technique  is  an  iterative  one  which  converges  after 
a  large  number  of  iterations,  but  which  has  tremendous  flexibility.  Ultima¬ 
tely,  it  should  be  possible  to  design  bulk  and  surface  wave  transducers, 
and  the  matching  network  as  one  package,  adjusting  all  parameters  for  the 
optimum  conditions,  whether  they  be  minimum  loss,  linear  phase  or  some 
other  parameter  like  good  transient  response,  which  is  of  interest  to  the 
designer.  The  technique  basically  minimizes  the  mean  squared  error  between 
a  desired  response  and  the  response  obtained;  so  it  has  a  great  deal  of 
flexibility. 

We  have  demonstrated  a  monolithic  Schottky  diode  devices  for  the 
first  time.  We  have  shown  that  v/e  can  indeed  make  Schottky  diode  cot  relators 
with  turn-on  times  of  the  order  of  a  few  ns  .  We  have  also  recently  been 
examining  the  possibility  of  putting  silicon  on  LiNb03  and  laser  an¬ 
nealing.  In  cooperation  with  Gibbons,  we  are  intending  to  work  toward  an 
alternative  Schottky  diode  device  based  on  the  silicon-lithium  niobate 
configuration.  We  have  new  ideas  for  putting  the  diodes  on  the  underside 
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of  the  silicon  and  we  will  be  testing  these  systems  along  with  ZnO  on 
Schottky  diode  correlators  In  the  next  few  months. 

We  have  had  some  difficulties  with  the  ZnO  deposition  because  of 
problems  with  the  Au  layers.  Unfortunately  there  was  a  change  in  the 
substrate  holder  and  the  method  of  mounting  the  substrate  in  the  holder, 
which  had  not  been  checked  out  sufficiently  carefully.  Consequently,  the 
thermal  contact  was  poor  and  the  gold  layers  that  were  being  put  down, 
which  had  been  reproducible  for  several  years,  v/ere  now  no  longer  properly 
oriented,  and  therefore  led  to  poor  ZnO  .  Consequently,  most  of  our 
transducers  had  rather  poor  performance,  even  though  the  quality  of  the 
ZnO  on  the  SiO£  was  excellent,  and  better  than  v/e  had  ever  seen  before. 

Unfortunately,  our  checks  with  a  bulk  wave  transducer  were  not  reliable 
because  the  bulk  wave  transducers  worked  well  with  the  use  of  a  different 
substrate  holder.  This  difficulty  has  now  been  corrected  and  we  should 
not  have  too  many  problems  in  the  future. 

B.  Research  Program  Plan 

We  will  continue  to  test  ZnO  deposition  on  a  wide  range  of  materials. 

We  will  carry  out  some  laser  annealing  experiments.  We  will  construct  and 
test  broadband  ZnO  on  Si  correlators  and  test  out  further  ZnO  on 
Si  Schottky  diode  correlators. 

C.  Major  Accomplishments 

We  have  demonstrated  a  Schottky  diode  ZnO  on  Si  correlator  for  the 
first  time.  We  have  arrived  at  entirely  new  techniques  for  designing 
broadband  ZnO  on  Si  devices  with  a  demonstrated  40  MHz  bandwidth. 

We  expect  to  be  able  to  obtain  still  larger  bandwidths  (60-70  MHz).  Sezav/a 
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mode  devices  and  important  new  matching  network  design  theories  with  broad 
scale  applications  have  been  developed. 


D.  Problems  Encountered 

No  major  problems  have  been  encountered. 


E.  Fiscal  Status 

Total  amount  of  contract 
Expenditures  &  commitments  through  12/31/79 
Estimated  funds  required  to  complete  work 
Estimated  date  of  completion  of  work 


$504,302 

$417,115 

$87,187 

30  September  1980 


F.  Action  Required  by  ARPA/ONR 

No  action  required  at  the  present  time. 
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II.  TECHNICAL  PROGRESS  REPORT 


A.  Sezawa  Wave  Oevices  and  Theory 

In  the  last  report,  theoretical  values  for  phase  velocity  and  coupling 
coefficient  for  (001)  ZnO  on  (001)  cut  (100)  prop.  Si  were  presented.  We 
have  since  made  Sezawa  wave  delay  lines  and  convolvers  which  exhibited  the 
high  coupling  coefficient  Av/v  =  .028  which  we  theoretically  predicted 
for  (001)  Si.  We  have  also  obtained  experimental  and  theoretical  results 
for  (111)  Si.  The  details  of  these  results  are  presented  in  Appendix  I, 
which  is  a  preprint  of  a  paper  submitted  to  Applied  Physics  Letters.  Only 
the  highlights  of  this  paper  and  some  additional  results  will  be  presented 
in  this  section. 

The  goal  of  the  present  work  on  Sezawa  waves  was  to  demonstrate 
broadband  thin  film  ZnO/Si  devices,  and  this  was  achieved.  Delay  lines 
with  3  dB  bandwidths  of  31  MHz  with  23  dB  loss  and  42  MHz  with  35  dB 
loss  were  demonstrated.  This  compares  very  favorably  with  the  best 
results  for  ZnO/Si  Rayleigh  wave  devices:  18  MHz  bandwidth  with  22  dB 
loss*  and  22  MHz  bandwidth  with  30  dB  loss. 2 

Sezawa  wave  convolvers  were  also  fabricated  using  (001)  Si.  These 
devices  all  had  a  high  insertion  loss  due  to  a  combination  of  two  problems: 
poorly  oriented  ZnO  under  the  transducers  and  interaction  losses  with  the 
electrons  in  the  central  Si  region.  However,  due  to  the  high  coupling 
coefficient  in  the  center  region,  these  devices  had  excellent  internal 
convolution  efficiency  (-27  dBm).  As  a  result  of  the  high  coupling 
coefficient  in  the  center  region,  the  insertion  loss  in  these  devices 
could  be  varied  by  45  dB  by  adjusting  the  top  plate  bias.  This  is  much 
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higher  than  the  1  to  3  dB  variation  seen  In  Rayleigh  wave  ZnO/Si  devices, 3 
and  indicates  that  high  gain  Sezawa  wave  amplifiers  could  be  made. 

Theoretical  values  for  coupling  coefficient  and  phase  velocity  have 
been  obtained  for  the  Sezawa  mode  (second  order  Rayleigh  mode)  for  ZnO 
on  substrates  of  (001)  Si,  (111)  Si,  and  for  (001)  LiNb03  on  (001)  Si. 
Theoretical  values  have  also  been  obtained  for  the  third  order  Rayleigh 
mode  for  ZnO  on  (001)  Si.  We  have  concluded  that  the  use  of  (001)  Si  as 
a  substrate  is  preferred  for  three  reasons:  (1)  the  highest  coupling 
coefficient  is  obtained  with  this  substrate;  (2)  the  peak  in  the  coupling 
coefficient  occurs  at  a  lower  value  of  hk  for  this  substrate,  and  hence 
it  has  the  largest  value  of  series  capacitance,  which  makes  tuning  easier 
and  allows  broader  bandwidth;  and  (3)  since  hk  is  smaller  for  (001)  Si, 
this  means  thinner  films  can  be  used,  and  we  expect  fewer  bulk  wave 
problems.  The  use  of  ZnO  on  (111)  Si  is  the  next  best  choice,  and  because 
of  the  lower  coupling  coefficient  and  higher  scries  capacitance,  the 
electronic  Q  is  twice  as  high  as  for  (001)  Si. 

Optimum  Matching  Network 

A  question  consistently  encountered  in  all  of  our  acoustics  work  has 
been  to  determine  what  is  the  maximum  bandwidth  attainable  for  a  given 
loss  if  the  ideal  matching  network  were  used.  The  answer  to  this  question 
allows  us  to  determine  if  our  matching  networks  are  close  to  ideal,  and 
it  would  also  allow  us  to  compare  the  use  of  different  materials  and  wave 
modes. 

This  question  was  first  answered  by  Bode^  for  a  parallel  R-C 
combination,  and  later  by  Fano5  for  series  R-L  and  several  other  slightly 
more  complicated  situations.  Their  results  will  be  applied  hero,  although 
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we  expect  that  the  maximum  bandwidth  attainable  will  be  slightly  less 
since  the  radiation  resistance  varies  with  frequency,  and  since  we  are 
neglecting  pad  capacity.  The  neglect  of  pad  capacity  is  a  good  approximation 
for  Sezawa  wave  IDTs  and  L i NbO 3  IOTs,  but  a  poor  approximation  for  ZnO/Si 
Rayleigh  wave  devices. 

For  an  IDT,  the  minimum  loss  L  for  a  given  fractional  bandwidth  f 
for  an  ideal  coupling  network  is 

1 

L  '  i  .  e-Wfi)  (1) 

where  Q  =  l/«RaoCy  and  Rao  is  the  radiation  resistance  at  center 
frequency  and  Cy  is  the  transducer  finger  capacity.  This  relationship 
is  shown  in  Fig.  1  for  a  number  of  interesting  cases.  It  can  be  seen  that 
ZnO/Si  Sezawa  v/ave  IDTs  and  LiNb03  Rayleigh  wave  IDTs  are  essentially 
equivalent  while  conventional  Rayleigh  wave  ZnO/Si  devices  have  inherently 
much  higher  electronic  mismatch  loss. 

These  results  indicate  that  40  MHz  bandwidths  for  our  Sezawa  wave 
IDTs  are  achievable  with  much  lower  insertion  losses  than  we  have  observed. 

Consequently,  we  developed  an  optimization  algorithm  which  will  determine 
the  optimum  component  values  of  an  arbitrary  matching  network.  This  algo¬ 
rithm  is  described  in  the  next  section. 

Optimization  Algorithm 

This  algorithm  is  based  on  our  earlier  studies  of  adaptive  filters. 

It  finds  the  component  values  which  minimize  the  mean  square  error 
between  the  desired  frequency  response  and  the  actual  frequency  response 
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over  a  given  frequency  range.  The  number  of  components  and  the  configu¬ 
ration  of  components  must  be  specified  by  the  user. 


Consider  a  generalized  network  as  shown  in  Fig.  2  where  the  transducer 
frequency  response  is  X(w)  and  the  frequency  response  of  the  network 
and  transducer  is 


Y(o>)  =  H(oj)  X{o) 


(2) 


X  and  Y  can  be  voltages,  powers,  or  some  other  parameter  of  interest. 
If  the  desired  frequency  response  is  D(a>)  ,  then  the  error  is 

E(co)  *  D(w)  -  Y(u)  *  D  -  H(w)  X(a»)  (3) 


Suppose  that  H  is  a  function  of  N  adjustable  real  parameters  such  as 
capacitors,  inductors,  or  resistors.  These  parameters  are  positive  and 
real  and  will  be  labelled  3n  .  Then 


EE*  =  (D  -  HX) (D*  -  H*X*) 

and  the  change  in  EE*  due  to  the  changes  in  A3n  is 


AEE* 


Using  Eq.  (4),  this  becomes 


(4) 


(5) 


L3H 

A3  Re  E*X  —  (6) 

n  *n 

We  wish  to  minimize  the  mean  square  error  over  the  frequency  range  from 
W1  w2  •  This  will  occur  if  we  ensure  that  the  change  in  the  mean 

square  error  is  negative,  i.e.,  choose 


(7) 


A6n 


dH 

E*X - du 


ae 


n 


where  an  are  positive  constants  which  determine  the  adaptation  rate. 
The  total  change  in  the  mean  square  error  is  then 


mj  n  L  cjj  n 


Using  the  choice  given  in  Eq.  (7)  for  A$n  t  this  algorithm  has  been 
implemented  with  good  results.  However,  the  problem  with  this  algorithm 
is  that  an  depends  on  the  type  of  component,  and  it  is  orders  of  magni¬ 
tude  different  for  capacitors,  inductors,  and  resistors.  This  problem 
can  be  avoided  by  noting  that  the  maximum  an  for  convergence  is 


a 


max 


/  EE*du 


Re 


[/ 


E*X 


_9H_ 


du 


(9) 


So  if  we  choose  an  =  Aannax  where  A  is  independent  of  n  and  of  the 
order  of  .1  to  .01  ,  then 


/  EE*du 

A3  =  A - 57-, - 

n  Re  J  E*X  du 

This  algorithm  v/i 1 1  find  the  optimum  component  values  where  optimum  means 
the  least  mean  square  error.  We  have  chosen  to  use  transferred  power 
(equivalently  insertion  loss)  for  Y  and  have  found  that  if  the  desired 
response  D  is  set  at  0  dR  loss,  then  a  flat  frequency  response  is 
obtained.  This  is  shown  in  Fig.  3  for  a  3  finger  pair  Sezawa  wave  IDT. 

If  the  desired  loss  is  sot  lower,  for  example  at  -2  dR  loss,  then  ripples 
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In  the  passband  are  obtained  (Fig.  3).  Thus  we  should  be  able  to  obtain 
43  MHz  bandwidth  with  only  3  dB  round  trip  electronic  mismatch  loss. 

The  original  experimental  and  theoretical  results  obtained  with 
4  finger  pair  Sezav/a  wave  IDTs  are  shown  in  Fig.  4.  Also  shown  in  that 
figure  are  the  improved  results  obtained  with  the  use  of  the  optimization 
program. 

Unbalanced  IDT  Excitation  on  Layered  Media 

For  Rayleigh  wave  devices  using  ZnO/Si  (Fig.  5a),  the  film  thickness 
is  -.006  .  The  maximum  bandwidth  (minimum  electronic  Q  )  is  obtained 

when  the  transducers  are  driven  in  a  balanced  mode.  These  transducers 
v/ere  often  used  in  an  unbalanced  mode  where  one  side  of  the  IDT 
is  driven  and  the  other  side  grounded  when  initial  measurements  were 
being  made  and  Av/v  determined.  The  measurement  in  this  case  is  simpler, 
since  a  balanced  to  unbalanced  transformer  is  not  needed.  The  Kino-VJagers 
theory  could  still  be  used  to  find  the  relation  between  Ra  and  Av/v 
since 

Ra  =  (R,  }/4  (10) 

unbal  dbal 


CT  =  2CT  (11) 

unbal  bal 

in  the  limit  of  A/h  0  . 

The  film  thickness  for  Sezawa  waves  is  typically  .25A  and  Eqs. 

(10)  and  (11)  no  longer  hold  since  the  capacitance  between  adjacent 
fingers  can  no  longer  be  neglected  (see  Fig.  5).  A  strong  reason  for 


using  Sezawa  wave  IOTs  In  an  unbalanced  node  is  that  the  unbalanced 
resistance  Is  of  the  order  of  150ft  instead  of  600ft  and  Cj  is 
~.7pF  instead  of  ~.35pf  .  Thus  matching  to  these  transducers  is  much 
easier  in  the  unbalanced  case. 

A  theory  has  been  developed  and  is  briefly  described  below  which 
predicts  Ra  and  Cj  for  arbitrary  voltages  Vj  and  V2  being  applied 
to  the  two  sides  of  the  transducer.  Special  cases  of  this  theory  are 
unbalanced  excitation  where  one  is  either  grounded  or  floating,  and 
balanced  excitation.  Comparison  of  theoretical  and  experimental  results 
are  also  included. 

We  shall  first  determine  the  capacity  of  the  structure  shown  in  Fig. 
5  and  then  the  radiation  resistance.  It  is  assumed  that  the  charge  is 
different  on  the  two  fingers  of  a  finger  pair,  but  that  every  finger  pair 
is  equivalent.  A  linear  approximation  will  be  used: 


V1  =  sllQl  +  s12Q2 


(12a) 


v2  =  S2iQi  +  S22Q2 


(12b) 


where  Sn  is  independent  of  Qi  .  The  derivation  is  simpler  if  we 
assume  the  finger  widths  and  spacings  are  the  same  for  all  the  fingers. 

In  this  case,  =  S22  and  S12  -  S21  •  The  more  general  case  is 
easily  solved  using  the  same  techniques.  The  two  unknowns  S\\  and  S\2 
can  be  found  by  considering  two  special  cases  where  the  relationship 
between  Vj  and  Qj  can  be  explicitly  calculated. 

Balanced  excitation  is  the  first  special  case  considered: 


Q2  =  -Ql 

V1  85  (SU  -  s12)  Ql 


ft! 


This  case  was  worked  out  by  Kino  and  Wagers6  with  the  result 


r^P»'V' 


1  (2V.) 

r  w  ~7>  =  2(sn  ‘  si2} 

Ls  yl 


1 


4 

irw 


^-^/sin  (2m+l)  1 

o'  (2ri+1)  f  (2m+1)Eo  +  Cl 


coth  yh 


(13) 


s  u  ■  -  i  - *  U  'P 

where  Cs  is  the  capacity  per  finger  pair,  (2Vj)  is  the  total  voltage 
across  the  capacitor,  d  is  the  width  of  the  finger,  and  £  is  the 
repeat  distance  of  the  transducer. 

For  the  second  special  case,  consider  the  case  when  one  side  of  the 
transducer  is  disconnected,  i.e.  Q2  -  0  .  Then 

V1  =  S11Q1 

As  shown  in  ref.  6,  the  finger  capacity  can  be  expressed  in  the  form 

c  =  — si— 

s  a 

w  |  o<f>dX 

To  find  the  relation  between  0  and  $  ,  it  is  easier  to  work  in  the 
Fourier  transform  domain  and  solve  Laplace's  equation  in  the  three 
regions  (substrate,  film,  and  air),  and  match  boundary  conditions.  The 
derivation  is  similar  to  that  of  ref.  6  except  that  there  is  a  term  in 
°k  and  <f>fc  with  no  spatial  frequency  variation,  and  the  solution 
of  Laplace's  equation  is  different  in  this  case.  The  result  of  the 
calculation  is 


’ll 


czzwt 


ity\L  £lx) _ 1 

™  '  "t-  '  c0 + 


coth  yh 
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Using  Eq.  (13),  Sj2  can  be  found 


*12 


1  ®  1  /ri*  2l!Wd  \“ 

1  1  /sin  -r-\ _ 

iiw  2m  \  '  c„  t  r  * 

(11=1  l  0  D 


czz«l 


coth  y.h 
p  '2m 


l_^p  _ J_  /sin 

*w“  (2m  +  1)  \  / 


™  m=0  {2m  + 


(2m+l )nd 
l 


Eo  +  cp  coth'*<2m*Dh 


Note  that  S-j 2  <  0  as  expected,  and  that  these  series  converge  very 
rapidly  (as  1/n3). 

We  can  now  calculate  the  capacitance  of  the  IDT  in  an  unbalanced 
mode  with  one  side  floating 


N 


or  with  one  side  grounded 


The  capacitance  in  these  cases  is  plotted  in  Fig.  6  as  a  function  of 
film  thickness  for  3  different  repeat  distances.  The  finger  widths 
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and  spacings  are  ft/ 4  In  this  calculation.  For  a  film  thickness  less 
than  ,25A  ,  it  doesn't  matter  whether  the  opposite  finger  is  floating 
or  grounded. 

To  derive  the  expression  for  radiation  impedances,  we  use 
the  expression  for  the  amplitude  of  the  forward  propagating  wave 


a(L)  =  A(Qi  -  Q2) 

where 


Nkd  kft 


The  radiation  inpedance  is? 

Za(L)a(L)* 


where 


and 


RaOF^ 

4|q/ 


(14) 
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In  the  balanced  mode,  Q]  ■  -Q2  and  Za  ■  RaoF(w)  .  Just  as  Kino 
and  Wagers'*  obtained. 

When  one  side  Is  grounded,  we  use  the  weak  coupling  approximation, 
i.e.,  Q  is  determined  from  electrostatic  considerations 


=  0  ■ 


S12Q1  +  S11Q2 


Thus,  using  Eq.  (14) 


Z 


a 


Raor(o,) 

4 


In  the  limit  h  0  ,  for  d/i,  «  1  ,  Sj 2  +  0  and 


Ra0F(») 

z  =  -12 -  „  .  0) 

a  m 


as  expected.  In  the  limit  h  -►  00  ,  then  S]2/Sn  -»■  1  and 


Za  *  RaOFU)  (h/x  »  1) 


v/h i c h  is  the  result  obtained  by  Kino  and  Wagers.® 

For  the  case  of  one  side  floating,  Q2  =  0  and 

7  „  R»0Fl",) 
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The  ratio  of  the  radiation  resistance  In  the  unbalanced  case  to  the 
balanced  case  Is  shown  in  Fig.  7. 

The  experimental  results  for  radiation  resistance  and  capacitance 
are  also  shown  in  Fig.  6.  The  agreement  with  the  theoretical  values  is 
good.  Note  that  the  radiation  resistance  is  increased  by  25%  when  the 
opposite  fingers  are  grounded  for  the  case  of  16  pm  periodicity  trans¬ 
ducers.  The  conclusions  of  this  theory  for  Rayleigh  wave  and  Sezawa 
wave  devices  are  given  in  Table  I.  Note  that  the  results  of  this  theory 
arc  important  for  the  Sezawa  wave  case,  but  unimportant  for  our  Rayleigh 
v/ave  devices  since  the  film  thickness  is  significantly  smaller. 

B.  Schottky  Hi  ode  Devices 

We  have  demonstrated  the  feasability  of  a  monolithic  Schottky  diode 
storage  correlator.  The  operation  of  this  device  is  similar  to  that  of  the 
pn  diode  storage  correlator,  with  the  advantage  of  having  the  ideal  diode 
transient  response  characteristic  of  the  Schottky  diode.  We  have  per¬ 
formed  experiments  that  indicate  that  it  is  indeed  possible  to  store 
signals  in  the  Schottky  diode  correlator  within  a  single  short  pulse  of 
a  duration  of  approximately  one  half  rf  cycle.  Figures  8  and  9 
indicate  the  autocorrelation  of  a  square  pulse  approximately  4.5  Psec 
long,  and  2  volts  in  amplitude.  In  Fig.  8,  a  6  nsec  7  volt  pulse  was 
applied  to  the  top-plate  to  store  the  pulse,  and  the  pulse  was  read  out 
approximately  50  ps  later  with  a  5  psec  ,  .75  volt  rf  modulated 
pulse  on  the  top  plate.  In  Fig.  9  the  output  correlation  of  the  stored 


signal  with  the  readout  signal  Is  observed.  We  have  found  that  a  new 
containment  box  fabricated  specifically  for  these  devices  has  been 
extremely  effective  in  reducing  the  rf  feedthrough  signals  that  so 
often  plague  these  devices. 

We  have  also  carried  out  work  on  improving  the  efficiency  of  these 
devices.  Apart  from  the  work  done  on  optimizing  the  deposition 
conditions  for  ZnO  films,  mentioned  elsewhere  in  this  report,  we  have 
investigated  some  of  the  other  loss  mechanisms  in  the  storage  correla¬ 
tor.  We  have  found  that  the  gold  film,  used  for  metallization  of  the 
top  plates  of  these  devices,  tends  to  give  excessive  loss,  this  is 
thought  to  be  due  to  scattering  from  grain  boundaries  in  the  polycrys¬ 
talline  gold  films.  These  losses  vary  over  a  wide  range  depending  on 
the  precise  deposition  conditions  of  the  gold  films,  but  are  typically 
of  the  order  of  3  dB  per  centimeter  with  a  1000A  gold  film.  These 
losses  are  far  above  the  theoretical  values  that  would  be  attributed  to  the 
viscosity  of  the  gold  film.  We  have  found  that  aluminum  top  plates  offer  a 
great  improvement  in  this  respect,  reducing  the  propagation  loss  to  essen¬ 
tially  negligible  values.  Improvements  in  insertion  loss  from  30  to  24  dB 
have  been  observed  by  merely  switching  to  aluminum  top  plates. 

When  the  correlator  is  operated  in  the  short  pulse  storage  mode,  the 
limitation  on  the  storage  time  is  given  not  by  the  characteristic  of  the 
diode  itself,  but  by  the  series  resistance  that  the  diode  sees  between 
itself  and  the  voltage  source. 


It  is  Important  that  the  RC  time  constant  of  this  resistance  in 
series  with  the  capacitance  of  the  top  plate  be  kept  to  a  minimum.  For 
this  reason,  we  have  begun  investigating  the  possibility  of  fabricating 
devices  in  n-type  epitaxial  layers  grown  in  highly  doped  n+  substrates. 
Preliminary  tests  indicate  that  the  required  epi  thicknesses  and  dopant 
concentration  levels  are  within  the  capabilities  of  already  existing 
deposition  facilities  at  Stanford,  and  we  are  now  developing  depth 
profiling  routines  to  confirm  that  the  desired  doping  profiles  have  been 
achieved. 

We  are  constructing  a  further  series  of  Schottky  diode  devices 
v/ith  lower  leakage  current  at  the  diodes  with  a  more  careful  use  of 
processing.  Leakage  current  in  these  Schottky  diode  devices  has  been  a 
major  problem  as  compared  to  the  equivalent  air-gap  devices. 

We  are  also  investigating  an  alternative  system  for  Schottky  diode 
devices,  making  the  use  of  polysilicon  deposited  on  lithium  niobate 
(LiNbOg),  which  is  then  laser  annealed.  Two  sets  of  experiments  will 
be  carried  out  in  cooperation  with  Prof.  Gibbons.  The  first  is  the 
deposition  of  Si 0£  on  Li Nb03  to  protect  the  lithium  niobate,  and 
then  the  deposition  of  polysilicon  on  the  Si Og  .  These  experiments 
will  be  run  to  determine  whether  we  can  indeed  deposit  polysilicon  on 
the  lithium  niobate  without  cracking  it.  If  we  can  do  this.  Gibbons  is 
confident  that  wo  can  develop  annealed  silicon  with  a  very  low  surface 
density.  Thus,  if  we  are  successful  in  producing  silicon  on  lithium 
niobate,  we  should  have  excellent  device  quality,  provided  the  acoustic 
losses  are  not  too  great.  Thus,  if  they  succeed  in  producing  silicon  on 
lithium  niobate,  we  will  go  forward  ininediat.cly  to  put  transducers  on 
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the  same  substrate  and  test  acoustic  loss  and  the  performance  of  the 
device  as  a  convolver. 

The  second  set  of  experiments  arc  concerned  with  whether  we  can  put 
Schottky  diodes  on  the  underside  of  the  silicon  so  as  to  make  a  storage 
correlator.  The  basic  idea  will  be  to  deposit  platinum  layers  on  the 
silicon  dioxide,  then  deposit  the  polysilicon,  thus  forming  Schottky 
diodes  in  the  annealed  silicon.  We  will  try  this  first  with  a  silicon 
substrate  as  the  idea  has  merit  for  its  own  sake. 

Our  ultimate  intention  is  to  make  a  high-quality  Schottky  diode 
storage  correlator  which  has  a  fast  turn-on  time.  It  is  our  aim,  of 
course,  to  make  a  monolithic  correlator.  So  far,  ZnO  has  appeared  to 
be  the  best  way  to  carry  this  out,  but  since  there  is  a  long  time  period 
between  runs,  there  is  time  without  use  of  much  additional  manpower  to 
carry  out  these  investigations  on  lithium  niobate.  Our  feeling  is  that 
it  would  be  unwise  to  neglect  the  laser  annealing  approach,  for  it  might 
provide  major  breakthroughs  in  the  surface  wave  field. 


C.  ZnO  Technology  Development 

Using  our  new  planar  discharge  sputtering  ZnO  system,  we  are  con¬ 
tinuing  our  characterization  of  ZnO  films  sputtered  at  high  substrate 
temperatures  up  to  450°C  .  Our  best  results,  as  measured  by  the  tuned 
round  trip  insertion  loss  of  a  sapphire  delay  line,  were  obtained  at  a 
substrate  temperature  of  450°C  .  The  following  table  summarizes  these 


results . 


r 


^sub 

F0 

I.L. 

450°C 

300  MHz 

3  dB 

500  MHz 

5  dB 

1  GHz 

7  dB 

2  GHz 

11  dB 

We  found  that  the  same  insertion  loss  is  obtained  with  different 
material  backing  layers  such  as  Au,  Pt,  and  A1 ;  as  long  as  these  metals 
are  well  oriented.  We  have  had  problems  depositing  well-oriented  ZnO 
on  Au  for  our  surface  wave  work.  This  problem  was  due  to  the  poor 
orientation  of  the  gold  film.  This  poor  orientation  of  the  gold  was  due 
to  an  unintentional  change  in  one  of  the  deposition  parameters  due  to  a 
change  in  the  substrate  holder.  This  difficulty  has  been  corrected  and 
we  should  be  able  to  make  surface  wave  devices  without  any  problems. 

We  have  carried  out  a  calculation  on  the  sputtering  profile  from  a 
planar  magnetron  discharge  sputtering  system,  to  determine  the  target 
*  substrate  spacing  for  maximum  uniformity  of  sputtered  fibers.  We  assume 

f  that  the  sputtering  from  the  erosion  ring  follows  a  cosine  law  distribution 

(a  valid  assumption  for  high-energy  sputtering).  Fig.  10  is  a  plot  of  the 
normalized  thickness  versus  distance  from  the  target  on  axis  (r  =  0)  . 

•  r 

Fig.  11  is  a  plot  of  the  thickness  profiles  of  the  sputtered  film  versus 
radial  distance  away  from  the  center  of  the  target,  for  various 
distances  zq  away  from  the  target.  The  thickness  is  normalized  to 
its  value  at  r  =  0  .  For  optimum  uniformity  of  the  sputtered  film,  the 
target  substrate  spacing  should  be  equal  to  1.3  times  the  radius  of 
the  erosion  ring.  It  is  obvious  from  Fig.  10  that  by  decreasing  our 
present  substrate  target  spacing  from  7.5  cm  to  5  cm  ,  wo  could  double 

-19- 


our  sputtering  rate  and  improve  the  uniformity  to  better  than  1%  over 
a  region  with  a  diameter  of  6  cm  .  We  are  presently  setting  up  an 

In  order  to  obtain  the  optimum  sputtering  conditions,  we  have  started 
a  detailed  study  to  evaluate  the  quality  of  ZnO  thin  films  sputtered 
under  different  conditions,  using  x-ray  diffractometry,  R.E.D,  S.E.M  on 
the  surface  and  on  a  fracture  edge,  and  resistivity  measurements. 
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TABLE  I 


Summary  of  Transducer  Parameters  for  Present  Devices 


Type 

Mode 

R/(RaOF(«)J. 

-C/Cbal 

Sezawa  Wave  Devices 

balanced 

1 

1 

(h  =  8  pn  ,  d  =  32  pm  , 

floating 

.25 

1.75 

d  =  8  ym) 

grounded 

.31 

1.80 

Rayleigh  Wave  Devices 

balanced 

1 

1 

(h  =  8  pm  ,  d  =  32  pm  , 

floating 

.25 

2 

d  =  8  pm) 

grounded 

.25 

2 
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FIGURE  CAPTIONS 


1.  Minimum  loss  versus  bandwidth  for  an  ideal  tuning  network. 

2.  Generalized  network  and  load  (IDT)  used  in  least  mean  square  optimiza¬ 
tion  algorithm. 

3.  Insertion  loss  versus  frequency  for  a  3-finger  pair  Sezawa  wave  IDT  for 
a  network  which  gives  minimum  ripple  (solid)  and  a  network  which  gives 
maximum  3  dB  bandwidth. 

4.  Insertion  loss  versus  frequency  for  a  4-finger  pair  Sezawa  wave  IDT 
showing  the  original  theoretical  and  experimental  results,  and  re¬ 
sults  of  using  the  optimization  program. 

5.  Schematic  drawing  of  a)  device  with  h/X  =  .006  (Rayleigh  wave 
device);  and  b)  device  with  h/X  =  .25  (Sezawa  wave  device). 

6.  Capacitance  versus  film  thickness  for  unbalanced  transducers, 
a)  w  =  1  mm  ;  b)  w  -  6  mm  . 

7.  Ratio  of  the  radiation  resistance  in  the  unbalanced  and  balanced 
cases  versus  film  thickness. 

8.  Pulse  used  to  store  acoustic  signal  in  storage  correlator: 

10  nsec  per  horizontal  division 

2  volts  per  vertical  division. 

9.  Storage  correlation  output  with  4  dB  of  gain  in  system: 

.02  volts  per  vertical  division 

2  ysec  per  horizontal  division- 

10.  Normalized  thickness  versus  distance  from  target  on  axis  (r  =  0)  . 

11.  Thickness  profiles  of  sputtered  film  versus  radial  distance  away  from 
center  of  target. 
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BROADBAND,  EFFICIENT,  THIN  FILM  SEZAWA  WAVE  IDTS 
J.  E.  Bowers,  B.  T.  Khuri-Yakub,  and  G.  S.  Kino 


Abstract 

The  phase  velocities  and  coupling  coefficients  for  Sezawa  waves  in 
(001)  ZnO  on  (001)  cut,  (100)  prop.  Si  and  on  (111)  cut,  (112)  prop.  Si  are 
presented.  Experimental  results  from  Sezawa  wave  delay  lines  are  presented 
which  represent  for  the  first  time  bandwidths  in  excess  of  30  MHz  in  a 
thin  film  monolithic  device. 
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Acoustic  surface  wave  (ASW)  signal  processing  devices  such  as  con¬ 
volvers  and  storage  correlators  can  be  monolithic,  thin  film  devices 
(typically  ZnO  on  Si)l  or  hybrid  arrangements  where  a  piezoelectric 
substrate  (typically  LiNbOj)  is  pressed  against  a  semiconductor  (typically 
Si).2»3  The  monolithic  approach  has  the  advantages  of  ruggedness,  smaller 
size,  and  fewer  spurious  signals,  and  the  devices  are  fabricated  using  stan¬ 
dard  planar  processing  and  are  potentially  much  cheaper  to  make.  The 
very  important  limitation  of  monolithic  ASW  devices  has  been  the  inherent 
narrow  fractional  bandwidth  (10-15%)  v/hich  is  a  result  of  the  low  coupling 
coefficient  (av/v  .004)  attainable  with  thin  film  Rayleigh  wave 
devices.  Larger  first  order  Rayleigh  wave  coupling  coefficients 
(av/v  _<  .015)  are  attainable  with  thicker  films  (h  ~  .5*  rather 

than  h  ~  .007a). 

Armstrong  and  Crampin'0  have  presented  theoretical  calculations  for 
the  coupling  coefficient  of  second  order  Rayleigh  waves,  commonly  called 
Sezawa  waves6  for  (001)  ZnO  on  (001)  Si.  They  found  a  peak  coupling 
coefficient  of  .024  at  a  velocity  of  4921  m/s  .  Elliot  et  al.?>8 
have  made  Sezawa  wave  signal  processing  devices  with  13%  bandwidth. 

In  this  paper,  we  present  theoretical  and  experimental  values  for 
the  phase  velocity  and  coupling  coefficient  for  Sezawa  waves  in  (001) 

ZnO  on  (001)  and  (111)  Si.  Information  is  also  presented  on  the  third 
order  Rayleigh  mode  on  (001)  Si.  Finally,  results  are  presented  on  Seza wa 
wave  delay  lines  with  fractional  bandwidths  of  19  and  25%  at  a  center  fre¬ 
quency  of  165  MHz. 

The  dependence  of  phase  velocity  on  the  thickness  h  of  the  different 
surface  acoustic  waves  that  are  possible  in  a  film  of  (001)  ZnO  on  a  (001) 
cut,  100  prop.  Si  substrate  are  shown  in  Fig.  1.  The  Sezawa  mode  is 


-  1  - 


cut  off  at  the  silicon  shear  wave  velocity,  which  is  5800  m/s  for  this 
orientation  and  propagation  direction.  This  mode  is  progressively  more 
lossy  as  k  is  decreased  below  the  cut-off  value  (kc  =  .73)  .  The 
cut-off  value  varies  by  as  much  as  15%,  depending  on  the  electrical 
conditions  at  the  surface  and  the  film-substrate  interface. 

The  coupling  coefficients  of  the  Sezawa  mode  are  shown  in  Fig.  2  for 
the  four  possible  transducer  configurations.  Note  that  we  predict  coupl¬ 
ing  coefficients  as  large  as  .028  ,  which  is  six  times  larger  than  the 
first  Rayleigh  wave  peak  for  this  orientation.^ 

The  coupling  coefficient  Av/v  must  be  calculated  at  a  constant 
value  of  u  ,  not  k  as  in  Armstrong  and  Crampin®,  since  the  theories^ >9 
for  transducer  impedance  all  use  Av/v  at  constant  to  .  Our  calcula¬ 
tions  for  Av/v  at  constant  k  agree  with  the  results  of  Armstrong  and 
Crampin5  and  yield  a  value  for  the  peak  coupling  coefficient  which  is  17% 
lower  than  the  correct  value.  Armstrong  and  Crampin  report  a  phase 
velocity  of  4921  m/s  ,  which  is  significantly  lower  than  our  theoretical 
and  experimental  results  and  the  experimental  results  of  Elliot  et  al J 
The  first  order  Rayleigh  wave  coupling  coefficient  is  .0002  at  the 
optimum  film  thickness  for  Sezawa  wave  coupling  (hk  =  1.65)  for  trans¬ 
ducers  at  the  free  surface  and  a  metal  shorting  plane  at  the  interface. 

Thus,  the  spurious  Rayleigh  wave  generation  is  50  dB  smaller.  The 
coupling  coefficient  for  Love  waves  to  these  transducers  is  zero,  so 
those  waves  are  not  excited. 

The  maximum  coupling  coefficient  for  third  order  Rayleigh  waves  was 
found  to  be  .0072  at  hk  -  5.8  .  This  mode  is  cut  off  at  hk  =  3.2 
and  is  not  desirable  for  these  applications  because  of  the  low  coupling 
coefficient. 
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The  Sezawa  wave  phase  velocity  and  coupling  coefficient  for  (111) 
cut  (112)  prop.  Si  are  shown  in  Figs.  3  and  4,  respectively.  The 
silicon  shear  wave  velocity  is  lower  (5000  p/s)  for  this  orientation 
and  propagation,  and,  consequently,  the  cut-off  occurs  at  a  higher  value 
of  hk:  2.15  . 

l/e  have  built  Sezawa  wave  delay  lines  on  both  the  (001)  and  (111) 
orientations.  The  experimental  agreement  with  theoretical  values  for 
phase  velocity  (Figs.  1  and  3)  and  coupling  coefficient  (Fig.  2)  is 
good. 

Using  the  coupled  resonator  tuning  network  shown  in  Fig.  5,  a  3  dB 
bandwidth  of  31  MHz  with  22.8  dB  insertion  loss  was  achieved  (Fig.  5). 
The  device  parameters  are  given  in  Table  I.  The  theoretical  plot  of 
insertion  loss  was  obtained  using  the  theoretical  values  for  phase  velocity 
and  av/v  and  the  measured  values  given  in  Table  I  for  the  device  and 
network  parameters.  The  equations  for  transducer  impedance  and  capaci¬ 
tance  derived  by  Kino  and  Wagers*'  for  balanced  excitation  were  modified 
to  apply  to  the  case  of  unbalanced  excitation.  Using  a  2-stage  it 
network, bandv/idths  as  large  as  42  MHz  were  attained,  although  the 
insertion  loss  was  higher  (35  dB). 

In  conclusion,  theoretical  and  experimental  values  for  the  phase 
velocities  and  coupling  coefficients  are  presented  here.  Bandwidths  in 
excess  of  30  MHz  were  achieved  for  the  first  time  in  a  monolithic  ZnO 
on  Si  device. 
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TABLE  I 


Delay  line  Parameters 


Substrate  orientation: 

(001) 

Ll  =  .040  pm 

Number  of  finger  pairs: 

4 

L2  =  .048  pm 

Wavelength:  32  urn 

L3  =  .9  pm 

Film  thickness:  S  pm 

Cj  =  36  pF 

Beam  width:  1  mm 

C2  =  40  pF 

Pad  capacitance:  .2  pF 

C3  =  9.6  pF 

- XC  1  "  1'.  _W.>  J 


FIGURE  LEGENDS 


Fig.  1  Phase  velocity  versus  normalized  filter  thickness  of  (001)  ZnO 
on  (001)  cut  (100)  prop.  Si. 

Fig.  2  Sezawa  wave  coupling  coefficient  versus  normalized  film  thick¬ 
ness  for  (001)  ZnO  on  (111)  cut,  (112")  prop.  Si. 

Fig.  3  Phase  velocity  versus  normalized  film  thickness  for  (001)  ZnO 
on  (111)  cut  (112)  prop.  Si. 

Fig.  4  Seza wa  wave  coupling  coefficient  versus  normalized  film  thick¬ 
ness  for  (001)  ZnO  on  (111)  cut,  (112)  prop.  Si. 

Fig.  5  Insertion  loss  versus  frequency  of  a  Sezawa  wave  delay  line  on 
(001)  cut  (100)  prop.  Si. 
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APPENDIX  II 


ADAPTIVE  DECONVOLUTION  USING  AN  ASW  STORAGE  CORRECTOR 

J.  E.  Bowers,  G.  S.  Kino,  D.  Behar,  and  If.  Olaisen 

Ginzton  Laboratory 
Stanford  University 
Stanford,  California 

Abstract 

A  new  analog  adaptive  filter  for  deconvolving  distorted  signals  is 
described  in  this  paper.  The  filter  uses  a  storage  correlator  which 
implements  a  clipped  version  of  the  LMS  algorithm  and  uses  a  special 
iterative  technique  to  achieve  a  fast  convergence.  The  new  filter 
has  a  potential  bandwidth  of  100  MHz  and  would  eventually  handle  pulsed 
signals  of  5  ysec  width.  For  signals  with  time-bandwidth  product  of 
less  than  100,  the  adaptation  time  is  less  than  1  msec  which  allov.s 
operation  in  real  time  for  most  applications  including  resolution  of 
radar  signals  in  a  cluttered  environment,  removal  of  echoes  from  tele¬ 
vision  signals,  deconvolution  of  distorted  signals  in  nondestructive 
evaluation,  and  also  in  telephony.  The  filter  is  particularly  suited 
for  radar  and  communications  as  it  processes  signals  directly  in  the 
VHF  range. 

Two  experiments  related  to  ghost  suppression  of  a  pulse  and  to  the 
field  of  NDE  are  described  in  this  paper.  The  results  arc  in  good 
agreement  with  computer  simulations  and  show  a  ghost  suppression  of 
15  dll  for  the  first  example  and  a  sidclobc  suppression  of  8  dll  for  a 
transducer  signal.  The  adaptation  time  is  less  than  450  nsec. 
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I.  INTRODUCTION 


Adaptive  filtering*  is  useful  in  removing  distortion  from  signals, 

particularly  when  the  distortion  varies  in  time.  Adaptive  filters  have 

been  used  to  perform  deconvolution  of  a  distorted  echo  pulse  in  an 

2 

acoustic  imaging  system  ,  to  equalize  the  distortion  xn  a  telephone 

3  4 

channel  ,  and  to  suppress  an  interfering  signal  . 

Most  adaptive  filters  have  been  implemented  using  digital  techniques. 
The  limitation  of  the  digital  approach  is  the  limited  bandwidth  (typi¬ 
cally  5  MHz)  and  the  practical  limit  on  the  number  of  taps  dictated 
by  the  complexity  and  power  consumption^- 

An  analog-digital  hybrid  approach  has  been  implemented  using  KOS  LSI 
technology^.  This  has  the  advantage  of  lowering  the  power  consumption 
and  allowing  32  taps  to  be  used  without  undue  external  complexity. 

Large  dynamic  range  was  obtained  with  this  technique  (60  dB),  but  the 
bandwidth  was  limited  to  less  than  1  MHz. 

Most  analog  implementations  of  an  adaptive  filter  have  been  made 

7 

using  CCDs  with  analog  tap  weights  held  in  sample  and  hold  circuits  . 

The  limitations  here  are  the  narrow  bandwidth  and  variation  across  the 
chip  in  gain  and  threshold  levels.  The  alternative  approach  used  at 
Hughes  for  implementing  a  wideband  adapting  filter  is  to  employ  tapped 

ASW  filters  with  complex  computer-controlled  systems  for  adjusting  the 

*  •  ««.  8 
tap  weigr.ts  . 

We  shall  describe  here  a  relatively  simple  all-analog  approach  to 

adaptive  filtering  which  uses  the  least  mean  squares  (LM5)  algorithm 

to  find  the  optimal  set  of  tap  weights.  An  ASW  monolithic  ZnO/Si  storage 

g 

correlator  with  8  MHz  bandwidth  and  the  equivalent  of  24  taps  was  used  . 
The  advantages  of  this  approach  arc: 
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(1)  fast  iteration  rate  (100  kllz)  winch  means  a  short  learning 
time  (-100  ys)  and  good  ability  to  track  time  varying  signals; 

(2)  simple  external  connections; 

(3)  large  potential  bandwidth  (100  MHz); 

(4)  large  potential  number  of  taps  (1000); 

(5)  low  power  consumption  (.1  W) ; 

(6)  it  is  suited  for  radar  and  communications  systems  because  it 
can  operate  directly  at  the  IF  frequency. 

A  disadvantage  is  that,,  at  the  present  time,  the  dynamic  range  is  less 
than  for  digital  systems. 

A  brief  introduction  to  the  I.MS  algorithm  is  provided,  followed 
by  a  new  derivation  of  the  LMS  algorithm.  As  the  class  of  applications 
discussed  here  utilizes  continuous,  quasistatic  filters,  a  derivation 
in  the  frequency  domain  rather  than  the  time  domain  allows  more  physi¬ 
cal  insight  into  the  processing  capability  of  the  LMS  algorithm  and  into 
its  advantages  and  limitations  for  these  kinds  of 
applications. 

Two  adaptive  deconvolution  experiments  which  used  the  storage  corre¬ 
lator  arc  described.  In  the  first  experiment,  an  undesirable  time- 
delayed  ghost  pulse  is  removed.  In  the  second  experiment,  the  ringing 
in  the  pulse  response  of  an  acoustic  transducer  is  removed  by  adaptive 
filtering.  The  results  of  these  experiments  are  compared  to  computer 
simulations  with  the  LMS  algorithm.  Conclusions  arc  drawn  regarding 
which  properties  and  nonli.nearitics  of  the  storage  correlator  arc 
limiting  the  performance  of  t lie  adaptive  filter. 


II.  EXPERIMENTAL  PROCEDURE 


A.  The  l.MS  Algorithm 

Consider  an  adaptive  filter  with  an  input  x(t)  and  output 
y(t)  =  x  *  w  .  The  tap  weights  w  arc  adjusted  after  each  iteration 
such  that  the  output  y(t)  converges  to  a  desired  signal  d(t)  .  The 
time  it  takes  for  convergence  to  occur  is  commonly  called  the  learning 
or  training  time.  After  the  optimum  set  of  weights  w  is  determined, 
the  filter  can  be  used  to  remove  the  distortions  in  signals.  For 
example,  in  a  TV  system  with  a  ghost,  the  filter  can  be  trained  on  the 
sync  signal  (which  occurs  at  the  end  of  each  line)  to  remove  the  echo 
sync  signal,  and  then  the  entire  TV  line  can  be  passed  through  the 
filter,  and  the  "ghosts"  will  be  removed. 

The  error  e(t)  between  the  desired  signal  d(t)  and  the  output  of 
the  filter  y(t)  is 


e(t)  =  d(t)  -  y(t  ) 


(1) 


where 


y  =  x  *  w  (?.) 

and  *  indicates  convolution.  Then  the  LMS  algorithm  specifies  that 
to  minimice  the  error,  w  must  he  changed  by  Aw  such  that 

Aw  =  px  v-V  c  '[?>) 

where  -/V  indicates  correlation  and  p  ir.  a  constant  which  determines 
the  rate  of  convergence  of  the  system. 


Thus  it  follows  from  Eqs.  (1),  (2),  and  (3)  that 


Aw  =  yx  ■£' (d  -  x  *  w)  (4) 

A  proof  is  given  in  the  following  section  that  if  these  adjustments 
Aw  are  made,  then  the  weights  exponentially  approach  the  Wiener  filter 
solution;  this  is  a  filter  which  minimizes  the  average  mean  square  error, 
over  the  frequency  band  of  interest,  between  its  output  signal  and  a 
desired  signal  D(to)  . 

It  can  be  seen  from  the  expression  for  the  tap  weight  adjustment 
(Eq.  (4))  that  for  each  iteration,  a  convolution  with  x  ,  and  a  corre¬ 
lation  with  x  must  be  performed.  In  addition,  an  analog  storage  regis¬ 
ter  is  needed  for  w  ,  and  it  must  be  possible  to  adjust  each  of  the  tap 
weights. 

The  ASK  storage  correlator  is  ideally  suited  for  this  purpose,  as  it 
can  perform  all  of  these  operations  (convolution,  correlation,  addition, 
and  storage).  A  big  advantage  of  this  approach  is  its  simplicity.  All 
of  the  necessary  operations  are  performed  inside  a  two  (or  three)  port, 
device  in  real  time. 

B.  The  Storage  Correlator 

The  ASK  storage  correlator  is  a  device  which  ullizes  the  interaction 
between  a  propagating  acoustic  surface  wave  and  a  charge  distribution 
stored  in  an  array  of  diodes  to  store,  correlate,  or  convolve  two  broad 
bandwidth  signals  in  real  time.  A  series  of  signals  may  be  iddcd 
together  by  successively  storing  them  in  the  diode  array.  A  Irief  descrip¬ 
tion  of  the  storage  correlator  follows.  Complete  descriptions  of  the 
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device  and  its  operation  os  well  as  a  theory  of  the  device 

are  given  in  the  literature. 


A  schematic  drawing  of  a  storatc  correlator  is  shown  in  Fig.  1.  If  a 
modulated  carrier  x  is  applied  to  the  acoustic  port  (#2)  and  a  signal  e 
to  port  #3  ,  then  the  charge  distribution^ 

Aw  =  x  -fe  c  (5) 

is  added  to  the  charge  distribution  w  stored  in  the  diode  array.  If 

a  signal  x  is  then  applied  to  port  #3,  then  the  convolution 

y  =  x  *  w  (6) 

% 

is  the  output  at  port  # 2.  The  LMS  algorithm  can  be  realized  by  repeating 
the  process  many  times.  Eventually,  the  tap  weights  w  converge  to  an 
approximation  to  the  Wiener  solution,  and  the  programming  is  complete. 

The  filter  may  then  be  used  to  deconvolve  other  signals  which  have  been 
distorted  in  the  same  manner. 

The  external  cop’  actions  used  to  operate  the  storage  correlator 
arc  indicated  in  Fig.  2a.  A  much  simpler  method  is  indicated  in  Fig.  2b 
when  two  acoustic  ports  of  the  correlator  are  used,  and  the  external 
delay  line  and  the  switches  are  not  needed.  A  difference  amplifier. is  not 
needed  to  subtract  d(t)  and  y(t.)  ;  rather  the  carrier  frequency  is 
adjusted  so  that  there  was  a  phase  difference  of  it  between  the  delayed 
signal  y(t)  and  the  desired  signal  d(t)  .  This  is  explained  more 
fully  in  the  next  section. 

The  storage  correlator  is  not  an  exact  implementation  of  the  Widrow 
type  LMS  adaptive  filter  for  two  reasons.  First,  the  LMS  algorithm 
requires  that  each  tap  he  changed  individually  and  the  error  signal  for 


the  next  tap  be  calculated  after  the  previous  tap  is  changed.  .  However, 
in  our  implementation,  the  entire  error  signal  is  calculated,  and  the 
entire  set  of  taps  is  changed  at  once.  This  difference  is  not  signifi¬ 
cant  except  for  the  case  of  very  high  gain  and  very  fast  convergence. 

The  second  difference  is  that  the  effect  of  the  plate  signal  (port  #3) 
on  the  tap  weights  is  not  linear  (Tig.  3).  If  the  error  signal  is  less 
than  a  threshold  determined  by  the  characteristics  of  the  devices  (4  V), 
the  tap  weights  are  unaffected.  The  error  signal  is  clipped  around  9V, 
so  only  ~6  dB  of  dynamic  range  is  available  at  port  U3.  If  the  error 
signal  is  fed  into  an  acoustic  port  (as  in  Fig.  4b),  then  at  least  35  dB 
dynamic  range  should  be  available. 

The  weight  signal  w(z)  is  stored  in  2000  diodes;  however,  the 
number  of  equivalent  taps  is  much  less  since  the  bandwidth  is  8  MHz 
and  t lie  maximum  signal  duration  is  3  ps  .  Consequently,  the  time  band¬ 
width  product  is  24.  The  device  was  represented  in  the  computer 
simulation  by  a  24  tap  transversal  filter.  The  necessary  correlations 
and  convolutions  were  calculated  in  the  computer,  and  the  error  signal 
was  modified  according  to  the  transfer  characteristic  shown  in  Fig.  3. 

It  can  be  seen  that  this  transfer  characteristic  closely  approximates 
the  experimentally  observed  characteristic. 

III.  THEORY 
h. _ Introduction 

The  LMS  algorithm  has  been  analyzed  in  detail  by  Widrow^’*^.  The 
first  task  in  his  analysis  is  to  transform  from  the  coordinate  system 
consisting  of  N  samples  from  the  waveform  to  a  set  of  normal  coordi¬ 
nates  where  the  correlation  matrix  is  diagonal;  the  eigenvectors  of 
the  correlation  matrix  arc  mutually  orthogonal,  and,  most  importantly. 


the  decay  modes  arc  uncoupled.  Once  this  transformation  is  accomplished, 
expressions  for  the  decay  constants  ,  misadjustment ,  etc.  may  be 
derived. 

In  the  analysis  presented  here,  the  transformation  is  very  easily 
accomplished  by  taking  the  Fourier  transform  of  all  quantities.  Corre¬ 
lations  are  then  simple  products,  and  the  convergence  of  a  quantity  at 
a  given  frequency  is  independent  of  other  frequencies,  i.e.,  the  inodes 
are  uncoupled. 


B.  Convergence 

Consider  the  situation  inf  Fig.  2,  where  all  signals  are  rf  modulated 
and  the  carrier  frequency  is  tog  .  The  desired  signal  d  (t)  and  the  dis 
torted  signal  x  (t)  are  repetitive,  and  as  the  number  of  iterations  n 
increases,  the  only  change  in  these  signals  is  a  phase  change 

jw  nT 

dn(t)  =  d(t)  e  U  (7) 


where  T  is  the  iteration  time.  In  the  frequency  domain 


Dn(w)  =  D(o>)  c 


*VT 


(8) 


Similar  expressions  hold  for  the  distorted  signal  X(co) 
The  filter  output  is 


Yn(w)  -  X(u)Wn(u)  c 


*  V 


(!>) 


and  tin;  error  signal  is 
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En(w)  =  Dn(w)  +  Yn(w)  e 


-J“otd 


j6>0nT  j”q 

B  D(co)  c  +  X(aj)Wn(oj)  c 


j(on(nT  +  Tq  -  Td) 

(10) 


where  is  the  length  of  the  external  delay  line,  and  Tq  is  the 

time  between  the  appearance  of  x(t)  at  the  acoustic  and  plate  ports. 
The  weight  adjustment  is 


w(t)  =  2paX(t)  en(t) 


(U) 


Equivalently,  we  can  write 

h’n+i  (w)  =  wn(“)  +  2yEji(w)X*;(w) 

Wn+1(«)  =  AWn(w)  +  B  (12) 


where 


A  = 


*  *VT0 

1  +  2yX(w)X  («)  e  u 


V 


B  =  2yD(co)X*  (co) 


(13) 


Equation  (12)  lias  the  solution 


V  (co)  «*  (W.  -  W  )An  +  W 

11  v  1  W  W 


(14) 


where  IL  is  the  initial  weight  distribution  and 
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w 

w 


D(o.)  jo»0(T0 

s  - e 

X(w) 


(IS) 


Equation  (14)  converges  only  if  ] A |  <  1  .  The  fastest  convergence 
is  obtained  for  the  carrier  frequencies  where  the  error  signal  is  the 
subtraction  of  D  and  Y  ,  i.c.  when 


IT 


w. 


I  td  -  To  I 


(2p  +  1) 


(16) 


where  p  is  a  large  integer.  If  the  feedback  gain  y  is  sufficiently 
small 


1 

V  <  -  (17) 

XX* 

hen  the  convergence  requirement  ]a[  <  1  is  satisfied  for  carrier 
frequencies  in  bands  of  width  Aw  around  the  frequencies  given  by 
Equation  (16) ,  where 

2 

Aw  =  - arccos  (yXX*)  (18) 

IT1)  -  Tol 

It  can  be  seen  from  Equation  (17)  that  when  y  is  near  its  tipper  limit, 
the  carrier  frequency  must  be  close  to  one  of  the  values  given  by 
Equation  (16). 

If  the  feedback  gain  is  significantly  less  than  its  maximum  value,  ■ 
it  follows  that 


10  - 


A  «  1  ♦  2yXX*  e 


*VT0  ‘  V 


a  e 


-1/T 


I 

(19) 


where 


1/t 


2yXX* 


^VTo 


V 


(20) 


The  convergence  of  the  weight  vector  (Eq.  (14))  is  approximately  exponen¬ 
tial.  So  we  can  write 


W 


~  (W1  -  IV  ) 

1  w' 


-n/r 
e  + 


W 

w 


(21) 


The  error  is 


E  «  D  -  X*W 
n  n 

En  *  x*  O'-'i  ~  \)  c"n/T  (22) 

The  error  decays  exponentially  to  zero.  In  practice,  the  bandwidth  of 
the  system  is  limited  by  the  bandwidth  B  of  the  acoustic  transducers. 
Consequently,  the  error  outside  this  frequency  range  is  not  affected  by 
the  adaptation  process.  It  will  also  be  noted  that  phase  distortion  in 
the  system  is  cancelled  out  by  the  application  of  the  convolution  and 
correlation  process  in  turn.  So,  although  this  is  basically  a  feedback 
system,  there  is  no  problem  with  instability  provided  the  convergence 
criteria  arc  satisfied. 

h’c  have  implicitly  assumed  that  the  signals  x(t)  and  d(t)  are 
sufficiently  short  compared  to  the  time  length  of  the  diode  array  that 
all  of  the  signal  Afa(t)  is  stored  in  the  diode  array.  If  this  is  not 
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the  case,  then  different  frequencies  do  not  decay  independently,  and 
the  analysis  is  more  difficult.  An  important  result  of  this  more 
difficult  case  is  that  the  error  does  not  decay  to  zero  in  the  passband. 

If  noise  on  the  input  signal  X  is  included  in  this  analysis,  and 
if  the  noise  on  the  plate  and  acoustic  ports  arc  correlated,  as  they 
would  be  if  a  delay  line  were  used  to  generate  the  second  X  signal 
which  is  needed  during  each  iteration,  then  the  expression  for  the  con¬ 
vergence  of  the  weight  signal  is  again  given  by  Equation  (22)  where 

DX* 

W  =  -  (23) 

w  XX*  +  <NN*> 

1 

T  - - (24) 

2y(XX*  +  <NN*>) 

Equation  (23)  is,  of  course,  the  Wiener  filter  solution. 

IV.  RESULTS  AND  DISCUSSION 
A.  Echo  Suppression 

In  this  experiment,  a  square  pulse  0.4  ps  long  is  followed  by  an 
echo  pulse.  The  desired  signal  is  a  single  pulse  0.4  ps  long,  for 
an  echo  which,  is  6  dB  less  in  .amplitude  than  the  main  pulse  (I;ig.  4a), 
the  si.delobe  suppression  after  10  iterations  (200  ps)  is  15  dB,  as 
shown  in  Fig.  4b.  The  dependence  of  sidclobo  suppression  on  echo  height 
is  shown  in  Fig.  5.  The  results  of  computer  simulations  of  1  lie  LMS 
algorithm  with  22  taps  and  using  clipping  and  a  threshold  6  dB  below 
1  he  clipping  level  are  also  shown  in  I-'ig.  5.  The  computer  simulation 
agrees  very  well  with  experimental  results  except  that  the  maximum  side- 
lobe  suppression  is  23.5  dB  which  is  4.5  dB  higher  than  was  experimentally 
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obtained.  The  reason  is  that  the  spurious  signals  generated  in  the  device 
during  read  out  limits  the  dynamic  range  of  the  signal  input. 

An  important  result  obtained  in  this  experiment  is  that  spurious 
acoustic  signals  generated  by  the  plate  readout  signal  can  be  suppressed 
by  up  to  13  dB  as  a  result  of  the  adapting  process.  The  filter  does 
not  distinguish  between  echoes  and  distortions  generated  externally  or 
by  the  device  itself.  This  result  is  demonstrated  in  Fig.  6.  The  upper 
trace  is  the  adaptive  filter  result  after  removing  an  echo  from  a  0.4  ps 
long  pulse.  If  all  signals  except  the  plate  readout  signal  are  removed, 
the  filter  output  (lower  trace)  is  the  spurious  signal  generated  by  the 
plate  signal,  A  large  spurious  signal  can  be  seen  when  previously  there 
had  been  a  null. 

The  advantage  of  computer  simulation  of  this  adaptive  filter  is  that 
the  threshold  and  clipping  levels  may  be  easily  changed  to  see  what  effect 
they  have  on  the  performance  of  the  adaptive  filter.  These  results  are 
summarized  in  Table  I.  The  computer  simulated  linear  LMS  result  is 
given  in  the  first  row.  If  a  threshold  level  is  included,  then  the 
computer  simulation  converges  much  faster,  but  to  worse  results.  The 
values  obtained  arc  in  agreement  with  experimental  results.  If  the 
feedback  gain  is  increased  so  that  the  error  signal  is  now  clipped, 
then  much  better  sidelobe  suppression  is  obtained  in  both  the  experimental 
and  computer  simulation  cases.  We  note  that  clipping  of  the  signal 
increases  the  rate  of  convergence  radically,  as  has  been  noted  by  others. 
The  algorithm  employed  is  therefore  known  as  the  clipped  LMS  algorithm. 

Regardless  of  the  shape  of  the  desired  signal  (d)  and  input  signal 
(x),  it  was  always  experimentally  observed  that  the  feedback  gain  must  be 
large  enough  to  strongly  clip  the  error  signal  during  the  first  few 
iterations  for  the  optimum  filtering.  After  many  iterations,  the  error 
signal  is  only  slightly  clipped. 
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Tho  threshold  and  clipping  levels  can  be  changed  by  modifying  the 
design  of  the  storage  correlator.  For  this  reason,  a  series  of  computer 
simulations  was  made  to  determine  the  effect  of  having  other  threshold 


to  clipping  ratios  (0)  than  the  value  o  =  .5  which  exists  in  the 
present  devices.  Only  the  ratio  a  is  important  here.  The  minimum 
r.m.s.  error  does  not  change  if  the  threshold,  clipping,  and  gain  values 
are  all  multiplied  by  the  same  value.  The  asymptotic  value  of  the  r.m.s. 
error  is  shown  as  a  function  of  o  in  Fig.  7  for  the  case  of  a  6  dB 
echo.  The  optimum  value  for  the  gain  is  used  at  each  point.  The  end 
point  o  =  1.0  corresponds  to  the  situation  of  hard  clipping  where  the 
threshold  level  equals  the  clipping  level,  and  the  feedback  error  is 
cither  1  or  0.  The  r.m.s.  error  increases  linearly  with  o  .  This  is 
intuitively  expected  since  any  error  values  less  than  the  threshold 
do  not  affect  the  tap  weights  and  arc  not  adapted  to  zero.  The  r.m.s. 
error  would  then  linearly  increase  with  threshold  level.  Consequently,  the 
design  of  the  storage  correlator  should  minimize  the  threshold  value. 

B.  Reduction  of  Sulk  Transducer  Ringing 

The  object  of  this  experiment  was  to  improve  the  impulse  response  of 
an  acoustic  bulk  wave  1.25  Ml lz  transducer  used  for  acoustic  nondestructive 
evaluation.  The  impulse  response  is  shown  in  Fig.  8a.  The  desired  signal 
is  a  unipolar  pulse  with  a  width  of  .  3  ps  .  If  the  impulse  response  is 
correlated  with  itself  (Fig.  8b),  the  highest  sidclobc  is  1.7  dB  smaller 
than  the  peak.  However,  if  the  correlator  is  used  as  an  adaptive  filter, 
then  the*  sidclobc  level  can  be  reduced  to  7  dB  below  the  peak  after  10 
iterations  (Fig.  8c)  and  10  dB  below  the  peak  after  35  iterations. 

Tlie  growth  of  1  lie  peak  as  the  number  of  iterations  increases  is 
shown  in  Fig.  9.  The  computer  simulation  is  also  shown  and  can  he  seen 
to  be  in  good  agreement  with  the  experimental  result.  The  computer 
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simulation  is  in  good  agreement  with  the  experimental  result  except 
that  the  maximum  predicted  suppression  is  2.6  dB  better  than  was  experi¬ 
mentally  observed. 

Linear  LMS  theory  predicts*  that  the  time  averaged  mean  square 
error  (MSE)  should  decay  exponentially  to  a  constant  value.  The  ratio 
of  this  constant  value  to  the  Wiener  solution  is  called  the  misadjustment 
M  .  The  misadjustment  is  approximately  related  to  the  decay  constant  T 
by  the  relation* 


M  =  N/4t  (2s) 

where  N  is  the  number  of  tap  weights.  In  our  case,  N  is  the  time 
bandwidth  product  of  the  storage  correlator.  Thus,  the  faster  the 
convergence,  the  larger  the  misadjustment. 

The  experimentally  observed  decay  of  the  MSE  for  the  bull:  transducer 
case  is  shown  in  Pig.  9.  The  decay  is  exponential  except  for  the  first- 
few  iterations  when  the  error  signal  is  strongly  clipped  and  the  mean 
square  error  decays  faster  than  an  exponential.  The  experimental  decay 
constant  is  7.S,  and  the  misadjustment  calculated  from  Equation  (25)  is 
equal  to  .8  .  We  interpret  this  to  mean  that  the  final  result  is  close 
to  the  Wiener  solution. 

C.  Limitations 

The  bandwidth  of  this  device  (S  Mils)  is  larger  than  the  bandwidth  of 
other  methods  of  implementing  the  LMS  algorithm;  however,  larger  bandwidth 
are  often  desirable.  Our  technique  can  work  well  as  an  adaptive  filter 
for  a  1.26  MU?,  transducer,  but  it  works  poorly  for  a  2.5  Ml  lx  transducer. 
The  diode  array  length  (3  ]is)  is  also  a  limitation  since  we  cannot: 
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adapt  signals  which  arc  longer  than  this.  Also,  since  the  correlation 
signal  is  truncated,  the  filter  output  is  distorted.  These  two  limi¬ 
tations  are  manifestations  of  the  fact  that  if  the  correlator  had  a 
larger  time  bandwidth  product  than  24,  then  it  would  equivalently  have 
more  taps  and  could  adapt  a  broader  class  of  signals.  Time  bandwidth 
products  of  500  or  more  arc  possible  with  the  use  of  present  technology. 
Note,  however,  that  liquation  (25)  predicts  that  if  the  number  of  taps 
is  increased,  then  the  convergence  time  also  increases  (for  a  given 
level  of  misadjustment) .  Thus,  a  larger  number  of  taps  is  not  desirable 
for  applications  in  which  the  distortion  is  changing  rapidly.  For 
example,  if  the  timing  or  amplitude  of  the  echo  significantly  changes 
over  a  time  period  cf  200  ps  (10  iterations),  then  a  larger  number  of 
taps  would  not  be  desirable. 

V.  CONCLUSIONS 

It  was  demonstrated  that  a  -6  dB  echo  could  be  suppresed  to  -15  dB 
in  2  00  ps  (10  iterations).  The  sidelobo.  of  the  impulse  response  of  a 
1.25  MHz  acoustic  bulk  transducer  can  be  reduced  from  -2  dB  to  -10  dB 
with  adaptive  filtering. 

The  experimental  results  are  on  the  whole  in  excellent  agreement  with 
the  computer  simulations.  Thus  it  appears  that  the  storage  correlator 
is  being  operated  in  the  optimum  manner  and  that  the  only  major  distor¬ 
tion  introduced  by  the  storage  correlator  is  the  presence  of  the  threshold 
and  clipping  levels. 

There  arc  many  advantages  of  using  a  storage  correlator  to  implement 
the  clipped  1.MS  algorithm  for  adaptive  filtering.  Only  two  connections 
arc  needed  for  a  device  which  can  correlate,  convolve,  store,  and  add 
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broadband  rf  signals  in  real  time.  The  device  is  small  (1  cm  )  and 
has  the  potential  of  much  broader  bandwidths  (100  MHz)  than  is  obtain¬ 
able  with  digital  techniques.  The  device  is  far  faster  in  operation 
than  any  competitive  inverse  filter  with  its  bandwidth. 
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Experimental  Results  and  Computer  Simulations  for  a 
Filter  Input  of  a  .4  ps  Wide  Pulse  with  a  -6  dll 
Echo.  The  Computer  Simulation  Uses  24  Taps. 
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1.  Schematic  drawing  of  storage  correlator. 

2.  (a)  Actual  and  (b)  proposed  implementation  of  SAW  adaptive  filter. 

3.  Correlator  output  as  a  function  of  plate  (readout)  pulse  height. 

4.  (a)  Input  signal  and  (b)  output  signal  after  22  iterations. 

5.  Experimental  and  computer  simulation  results  for  echo  suppression 
experiment . 

6.  Echo  reduction  after  10  iterations  (upper  trace).  Spurious  signals 
generated  by  plate  read-out  signal  (lower  trace). 

7.  Minimum  r.m.s.  error  as  a  function  of  threshold  level. 

8.  (a)  Impulse  response  of  1.25  MHz  bulk  acoustic  transducer;  (b)  auto¬ 

correlation  of  impulse  response;  and  (c.)  filter  output  after  10 
iterations  of  adapting. 

9.  Output  pulse  level  during  adapting  process. 


0.  Mean  Square  Error  during  adapting  process. 
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